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The fluctuation-dissipation theorem (FDT) plays a fundamental role in understanding quantum 
many-body problems. However, its applicability is limited to equilibrium systems and it does in 
general not hold in nonequilibrium situations. This violation of the FDT is an important tool for 
studying nonequilibrium physics. In this paper we present results for the violation of the FDT in 
the Kondo model where the impurity spin is frozen for all negative times, and set free to relax 
at positive times. We derive exact analytical results at the Toulouse point, and results within a 
controlled approximation in the Kondo limit, which allow us to study the FDT violation on all time 
scales. A measure of the FDT violation is provided by the effective temperature, which shows initial 
heating effects after switching on the perturbation, and then exponential cooling to zero temperature 
as the Kondo system reaches equilibrium. 
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I. INTRODUCTION 

The fluctuation-dissipation theorem (FDT)i is of fun- 
damental importance for the theoretical understanding of 
many-body problems. It establishes a relation between 
the equilibrium properties of a system and its response 
to an external perturbation. In nonequilibrium situa- 
tions this powerful tool is in general not available: typical 
nonequilibrium situations are e.g. systems prepared in an 
excited state, or systems driven into an excited state by 
pumping energy into them. Since such nonequilibrium 
systems occur everywhere in nature, the investigation of 
nonequilibrium many-body physics has become one of 
the key challenges of modern theoretical physics. The vi- 
olation of the FDT in a nonequilibrium system plays an 
important role in such studies as it characterizes "how 
far" the system is driven out of equilibrium. 

Most widely investigated in this context are glassy sys- 
tems, that is systems with a very long relaxation time 
compared to the typical time scale of measurements. Due 
to the long relaxation times it is experimentally possible 
to measure the deviation from the FDT, i.e. to study the 
ageing effects: one observes a relaxation of the nonequi- 
librium initial stat e to wards equilibrium. For a review of 
this field see Refs. I2l3l 

However, these are classical systems at finite tempera- 
ture and therefore the classical limit of the FDT is stud- 
ied. Nonequilibrium zero temperature quantum systems 
provide a very different limit which has been studied very 
little in the literature (however, see Ref. 0). This is one 
of the main motivations for our work which looks at the 
FDT violation in the time-dependent Kondo model at 
zero temperature. Time-dependence is here introduced 
by freezing the impurity spin at negative times, and then 
allowing it to relax at positive times. Besides being of 
fundamental theoretical importance as the paradigm for 
strong-coupling impurity physics in condensed matter 
theory, Kondo physics is also experimentally realizable 
in quantum dots. The Kondo effect has been observed in 
quantum dot experiments 5 , and time-dependent switch- 



ing of the gate potential amounts to a realization of the 
time-dependent Kondo model^ which should be possible 
in future experiments. 

Our calculations here are based on recent work on the 
time-dependent Kondo model with exact analytical re- 
sults for the Toulouse point and results in a controlled ap- 
proximation in the experimentally relevant Kondo limit. 7 
We will see that the FDT is maximally violated at inter- 
mediate time scales of order the inverse Kondo temper- 
ature: the effective temperature becomes of order the 
Kondo temperature due to heating of the conduction 
band electrons by the formation of the Kondo singlet. 
The system then relaxes towards equilibrium and the 
FDT becomes fulfilled exponentially fast at larger times. 



II. FLUCTUATION-DISSIPATION THEOREM 

Consider an observable A which is coupled linearly to 
a time-dependent external field h(t). The Hamiltonian 
of the system is then given by 



H = Ho - h(t) A 



(1) 



where Hq is the unperturbed part of the Hamiltonian. 
The generalized susceptibility (or response) R{t,t') of the 
observable A at time t to the external small perturbation 
hit') at time t' is defined as 



R(t,t') 



S(A(t)) 



5h(t') 



(2) 



h=0 



Here (A(t)) = (A(t)) h -(A(t)) is the deviation of the ex- 
pectation value from its equilibrium value. If the system 
is in equilibrium before its perturbation by the field h, 
then R(t, t') depends only on the time difference r = t—t'. 
One introduces the Fourier transform of R(t) 



R(lj) 



R{T)e lu,T dT 



(3) 



2 



where the integration runs only over positive times as a 
consequence of causality. A simple calculation shows that 
the imaginary part of R{oj) is proportional to the energy 
dissipated by the system for a small periodic perturbation 
with frequency ui (see, for example, Ref. |8|) ■ Thus the 
response function determines the dissipation properties 
of the equilibrium system. 

For calculating the response function one defines the 
two-time correlation function 

C A , A (t, f ) = (A(t)A(t')) = ±Tr [A{t)A{t')p] , (4) 

with the operators in the Heisenberg picture 

Alt) = exp (iHt) A(0) exp (-iHt) . (5) 

The trace runs over all the states in the Hilbert space, 
p is the density matrix and Z the partition function. 
Symmetrized and antisymmetrized correlation functions 
Cs AA \ (t; t'), C\ AtA \(t,t') are defined in the same way 

C {A ,A}M = \({A(t),A(t')}) (6) 

C [A ,A]M = \{[A{t),A{t')}) . (7) 

The cumulant of the symmetrized correlation function is 

C[ c $ } (t,t') ee C {A , A} (t,t') - (A(t)) (A(f )) . (8) 

In the framework of linear response theory (that is for 
small perturbations) one then proves the famous Kubo 
formula— 



R(t,t') = 2i6(t-t')C [A>A] (t,t') 



(9) 



Here the time dependence of all operators is given by the 
unperturbed part of the Hamiltonian Hq. Since in equi- 
librium all correlation functions depend only on the time 
difference r = t — t', one defines their Fourier transform 
with respect to r 



r>( cum ) i, ,\ 



(10) 



If the initial state is the equilibrium state for a given 
temperature, i?(w) and cj™™j(w) are related by the fa- 
mous Callen-Welton relation 1 , which is also known as the 
Fluctuation-Dissipation theorem 



Iii.i = ta.nb (^) C{ C ™1M 



(11) 



Here (3 is the inverse temperature. For T — Eq. 
reads 



ImiJH^sgnHCf^jH. 



(12) 



Eqs. (jllfl and (|12fi » relate dissipation with equilibrium 
fluctuations, which is the fundamental content of the 
FDT. 



III. FDT VIOLATION IN NONEQUILIBRIUM 

Let us recapitulate why the FDT (II 1|) in general will 
not hold in quantum nonequilibrium systems. We will 
only consider the zero temperature case since it brings 
out the quantum effects most clearly; the generalization 
to nonzero temperatures is straightforward. 

We first consider how a typical experiment is actu- 
ally performed: the system in prepared in some initial 
state at time t = (not necessarily its ground state) and 
then evolves according to its Hamiltonian. A response 
measurement is then done by applying the external field 
after a waiting time t w > 0, and the response to this is 
measured a time difference r later. The Fourier trans- 
form with respect to (positive) time difference will then 
in general depend on the waiting time t w 



R(u,t w ) 



R(t w + T, t w )e^ T dr 



(13) 



Likewise in an experimental measurement of the corre- 
lation function the first measurement of the observable 
will be performed after the waiting time t w , and then at 
time t w + t the second measurement follows. From the 
experimental point of view this leads again to a one-sided 
Fourier transform 



(, , + 

{A, A} v ' " 



^(cum) / , 
^{A,A}\ 1 ^ 



t, t w ) cos(u)t)<1t 



. (14) 

If the system is prepared in its ground state, or if the 
system equilibrates into its ground state for sufficiently 
long waiting time ^ oo, then we can replace this one- 
sided Fourier transform by the symmetric version and 
arrive at the conventional equilibrium definition 1101) . 

However, for a nonequilibrium preparation at t = 
Eqs. (|13|) and (|14l) are the suitable starting point for the 
discussion of the FDT (|12|) . Let us therefore look at 
the FDT in the framework of lO and (|T1J- We follow 
the standard derivation- and introduce a complete set 
of eigenstates \n) of the Hamiltonian H, H\n) — E n \n). 
The matrix elements of the operator A are denoted by 
A nm = (n\A\m) in this basis. Then a matrix element of 
the susceptibility is given by 

R(t, t w )nn' = i0{t — tw) ^ A (15) 

m 

x ^ e i(E n -E m )t e i(E m -E n ,)t u _ e i(E n -E m )t w e i(E m -E n ,)t^ 

The imaginary part of l|13|) is 

lmR(uj, t w ) nn i — Re A 



(16) 



/o 



) 



with uj nm = E n — E m . For diagonal matrix elements 
ii = n' this implies 



Imi?(w, ty 



(17) 
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m 

Likewise for the correlation function 



(18) 



Jo 

and the diagonal matrix elements are 

{A A} ^ 1 w J nn — 



(19) 



If we take \n) = |GS) as the ground state of our Hamil- 
tonian, i.e. the system is in equilibrium, then we know 
Wnm = Eqs — E rn ^ and u) mn — E m — Eqs > for 
all m. For positive uj therefore only the first terms in 
(|17J) and (|19|l contribute, and for negative u> the second 
terms contribute: this just proves the zero temperature 
FDT (jT2)l with its sgn(w)-coefficient. 

Now let us assume the nonequilibrium situation de- 
scribed above where the system is prepared in some ar- 
bitrary initial state |NE) at t = 0. One can expand |NE) 
in terms of the eigenstates \n) of the Hamiltonian 



|NE)=^c n | 



(20) 



We allow for anisotropic couplings Ji = (J±, J±, J||) and 
consider a linear dispersion relation ek — vf k. We have 
studied two nonequilibrium preparations in Ref.0: I) The 
impurity spin is frozen for time t < by a large magnetic 
field term h(t)S z that is switched off at t — 0: h(t) ^> Tk 
for t < and hit) — for t > 0. II) The impurity spin is 
decoupled from the bath degrees of freedom for time t < 
(like in situation I we assume (S z (t < 0)) = +1/2) and 
then the coupling is switched on at t = 0: Ji(t) = for 
t < and Ji(t) — Ji > time-independent for t > 0. For 
both scenarios the impurity spin dynamics could then be 
described by an effective time-dependent resonant-level 
Hamiltonian in terms of fermionic solitons ^k consisting 
of spin-density excitations: 



E9kk>*l*A did - 1/2) ,*<Q 



H 



kk> 



, t > 
(23) 

with effective parameters gkk' and Vk (which also de- 
pend on scenario I or II). The impurity spin S z is given 
by S z — dJd — 1/2. Since the effective Hamiltonian 
is quadratic for both negative and positive times, it is 
straightforward to find an explicit solution for the im- 
purity orbital correlation functions and work out their 
dependence on the waiting time. Detailed inspection 7 
shows that the impurity spin dynamics is the same for 
both nonequilibrium initial preparations I and II. 



with suitable coefficients c„ . Then the relations IjRil) and 
(|18H are modified like 



Imi?(w,^) NE = 



m,n,n' 



m,n,n' 



where " . . ." stands for the same expressions as in <|16[) 
and (|18|l . In general this will lead to a nonzero difference 



sgn(w) C^^y(LJ,t w ) N E - Im R(uj,t w ) NE ^ ( 

and therefore the FDT is violated. We will next study 
the violation of the FDT explicitly in the time-dependent 
Kondo model, and in particular also show that the dif- 
ference (|21f) vanishes exponentially fast for large waiting 
times t w . 



21) 



V. TOULOUSE POINT 

At the Toulouse poin1*±I with J\\/2ttv f = 1 - 1/a/2 
the mapping to the effective resonant level model (|23|l is 
exact and the effective parameters gkk' and Vk are inde- 
pendent of k, k' . This allowed us to express the spin-spin 
correlation functions in closed analytical form 7 

^(cum) / , % def 

= \{{S z (t w ),S z {t w + t) }) - (S z (t w )) (S z (t w + r)) 



-2t/* e 



(!_ e -4WtB) 



(24) 



- (s{t) - s(t w + r) e -'»/* B + s(t w )e-^+^^y 
for the symmetrized part and 



IV. TIME-DEPENDENT KONDO MODEL 

We briefly review the results obtained in Ref. for 
the spin dynamics of the time-dependent Kondo model. 
The time-dependent Kondo model is described by the 
Hamiltonian 



\ " rSTJ 



a/3 

T " . 



k , a. 



a,0 



(22) 



C[Sz,Sz](T,t w ) d = ^ ( [S z (t w ),S z (t w +t)]) 



= -i. e - T/tB 



( s ( T )- s (^+ T ) e -^/*B 
+ S (^)e-^+^/ tB ) (25) 



del' 



for the antisymmetrized part. Here s(t) 
(^b/ 71 ") f Q °° dui sin(wr)/(l + w 2 tg) with the shorthand no- 
tation — TTwtK- w — 0.4128 is the Wilson number 



4 




FIG. 1: Universal curves for the spin-spin correlation func- 
tion Tk x C|g Um J y(uj,t w ) (solid line) and response function 
Tk x Im R(lu, t w ) (dashed line) at the Toulouse point. Notice 
the normalization of the equilibrium curve (t w — * oo) which 
follows from © with the operator identity S% = 1/4: this 
gives J °° (7r^g a i(w,t ro = oo)duj = tt/4. 



and tx = 1/Tk is the Kondo time scale, i.e. the inverse 
Kondo temperature. In this paper we use the defini- 
tion of the Kondo temperature Xk via the zero tempera- 
ture impurity contribution to the Sommerfeld coefficient, 

7imp = W1T 2 /3T K . 

From (|24|l and (|25|l one can obtain the Fourier trans- 
forms {T3J and ifHjl. Results for C c ^ m s Aw,t w ) and 
ImR(uj,t w ) for various waiting times t w are shown in 
Fig. 2] For zero waiting time t w = the FDT is triv- 
ially fulfilled since the system is prepared in an eigen- 
state of S z and therefore both functions vanish identi- 
cally, C{™ m ^ } (w,^ = 0) = lmR(w,t w = 0) = 0. For 
increasing waiting time the curves start to differ, which 
indicates the violation of the FDT in noncquilibrium. For 
large waiting time as compared to the Kondo time scale 
one can then see nicely that the curves coincide again, 
which shows that the system reaches equilibrium behav- 
ior for t w — > oo where the FDT is known to hold. From 
the curves in Fig. ^ one also notices that the maximum 
violation of the FDT occurs at zero frequency, while it be- 




0.001 



FIG. 2: Effective temperature T e s as a function of the waiting 
time t w at the Toulouse point. The inset shows the same 
curve on a linear scale to illustrate how fast the initial heating 
occurs. 



comes fulfilled more rapidly at higher frequencies. We in- 
terpret this as showing that high-energy excitations find 
"equilibrium-like" behavior faster than low-energy exci- 
tations probed by the small u> response. The high-energy 
components of the initial nonequilibrium state can "de- 
cay" more quickly for a given waiting time. 

At zero frequency Crg Un ^ Auj = 0,t w ) is non-zero for 
< t w < oo as if one were studying the spin dynamics of 
the equilibrium system at finite temperature. This leads 
to the definition of the effective temperature T e g via the 
zero frequency limit of 



lim 



Imi?(w) 



& = 0). (26) 



2T off 



This concept of an "effective temperature" is frequently 
used and well-established in the investigation of classi- 
cal nonequilibrium systems^ We suggest that it is also 
useful in a quantum noncquilibrium system by giving a 
measure for the "effective temperature" of our bath (i.e. 
the conduction band electrons) in the vicinity of the im- 
purity. 

We can see this explicitly by using H2(jfl to evaluate the 
effective temperature as a function of the waiting time; 
the results are shown in Fig. [5] One sees that the effec- 
tive temperature goes up very quickly as a function of the 
waiting time until it reaches a maximum of T c g m 0.4Tk 
at t w fts 0.1£k- After that the system cools down again. 
We can understand this by noticing that the conduction 
band is initially in its ground state with respect to the 
Hamiltonian for t < 0, therefore its effective temperature 
vanishes. As the spin dynamics is turned on at t = 
the Kondo singlet starts building up. Its nonzero bind- 
ing energy therefore initially "heats up" the conduction 
band electrons. After a sufficiently long time the Kondo 
singlet has been formed and then the process of energy 



5 



diffusion takes over: the binding energy that has initially 
been stored in the vicinity of the Kondo impurity diffuses 
away, the system equilibrates and the effective temper- 
ature goes back to zero. The behavior of the effective 
temperature therefore traces this competition of release 
of binding energy and energy diffusion away from the im- 
purity. Analytically one can show for very small waiting 
time t w <C £k 



1 



1 



w ln(7TUl t^/t w ) 



(27) 



and an exponential decay to zero temperature for long 
waiting time t w 3> £k 



T, 



_£ K e -W* K 



(28) 



Finally we want to emphasize that while the effective 
temperature seems a useful phenomenological concept 
for interpreting the u> = behavior in nonequilibrium, 
its definition l|26|l does not capture the small w-behavior. 
Since cj™"^ Au,t w ) ~ ^{s^s }( w = 0>*iu) ^ M is non- 
analytic for small u) and finite waiting time (see Fig. ^ 
and the discussion in Ref . @) , the long time decay of the 
spin-spin correlation function is always algebraic for all 
t w > (therefore characteristic of equilibrium zero tem- 
perature behavior), cj™"^ j(t, t w ) oc r~ 2 . 
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VI. KONDO LIMIT 

The Kondo limit with small coupling constants 
J±,J\\ — > is the relevant regime for experiments on 
quantum dots. In this regime the results in Ref. [7J are 
not exact, but were shown to be very accurate by com- 
parison with asymptotically exact results for t w = and 
r/tji ^ l- 13 For our purposes here the main difference 
from the Toulouse point analysis is the nontrivial struc- 
ture of the effective parameters gkk' and 14 from Ref. in 
(J2HJ. This makes it impossible to give closed analytic ex- 
pressions like (12411 or Q25JI. but the numerical solution of 
the quadratic Hamiltonian l|23|l is still straightforward. 
The results presented in this section were obtained by 
numerical diagonalization of l|23|l with 4000 band states. 
The numerical errors from the discretization are very 
small (less than 2% relative error in all curves)^ 

Fig. |21 shows the behavior of the spin-spin correlation 
function and the response function obtained in this man- 
ner. Similar to the Toulouse point results we observe 
a violation of the FDT for finite nonzero waiting time, 
< t w < oc. For t w — > oo one recovers the FDT ex- 
ponentially fast (|28|l as expected since the system equi- 
librates. The results for the effective temperature in the 
Kondo limit are depicted in Fig. 0] While the behav- 
ior of T e s(t w ) is somehow more complicated than at the 
Toulouse point, the interpretation regarding heating and 
cooling effects carries over without change. The main 
difference is that the maximum effective temperature is 



FIG. 3: Universal curves for the spin-spin correlation function 
Ik x C|g Um J y((v,tw) (solid line) and response function Tk x 
Im R(lu, t w ) (dashed line) in the Kondo limit (compare with 
Fig-EJ. 

already reached for t w w O.OMk in the Kondo limit. We 
interpret this as being due to the (dimensionful) bare cou- 
pling constants at higher energies that are larger than the 
renormalized low energy scale Tk and therefore lead to 
faster heating. 



VII. CONCLUSIONS 

Our investigation of the zero temperature quantum 
limit of the fluctuation-dissipation theorem in the time- 
dependent Kondo model provides some important lessons 
regarding its relevance in quantum nonequilibrium sys- 
tems. For the Kondo system prepared in an initial state 
with a frozen impurity spin, i.e. in a product state of sys- 
tem and environment, the FDT is violated for all nonzero 
waiting times t w of the first measurement after switch- 
ing on the spin dynamics at t = 0. For large waiting 
times as compared to the Kondo time scale the FDT 
becomes fulfilled exponentially fast, which indicates the 
quantum equilibration of the Kondo system. A quantita- 
tive measure for the violation of the FDT is provided by 
the effective temperature T e fji£ here defined via the spin 
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FIG. 4: Effective temperature T B g as a function of the waiting 
time t w in the Kondo limit. The line is a guide to the eye. 
The size of the datapoints (circles) indicates the numerical 
error. The datapoint for t w /tK = 5.5 is numerically indistin- 
guishable from zero. 

dynamics (|2bl and depicted in Figs.|21and0] It traces the 
buildup of fluctuations in the conduction band: Initially, 
the conduction band electrons are in equilibrium with 
respect to the Hamiltonian for t < 0. Then in the vicin- 
ity of the impurity they get locally "heated up" to T e g 
about 0.4T K (Toulouse point)/ 0.45T K (Kondo limit) due 
to the release of the binding energy when the Kondo sin- 
glet is being formed. Eventually, this excess energy dif- 
fuses away to infinity and T e g reaches zero again. In this 



sense the largest deviation from zero temperature equi- 
librium behavior occurs for t w rs O.Uk at the Toulouse 
point, and for t w w 0.03tK in the experimentally rel- 
evant Kondo limit, with very rapid initial heating (see 
the inset in Fig. [5}. These observations could be relevant 
for designing time-dependent (functional) nanostructures 
with time-dependent gate potentials^ since they give a 
quantitative insight into how long one needs to wait af- 
ter switching for the system to return to (effectively) zero 
temperature. 

From a theoretical point of view it would be interest- 
ing to study the FDT for other observables (like the cur- 
rent) and in other nonequilibrium quantum impurity sys- 
tems in order to see which of the above observations are 
generic. Notice that the "effective temperature" will gen- 
erally depend on the observable chosen for its definition 
in (|26t 15 . though we suggest that the qualitative behavior 
(rapid initial increase and exponential decrease) will be 
similar for all local observables. Work along such lines is 
in progress in order to substantiate the concept and no- 
tion of an "effective temperature" qualitatively charac- 
terizing the evolving nonequilibrium quantum state, and 
to explore its usefulness in quantum nonequilibrium mod- 
els in general. 
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